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Abstract 

We study non-linear contributions to the power spectrum of the curvature perturba- 
tion on super-horizon scales, produced during slow-roll inflation driven by a canonical 
single scalar field. We find that on large scales the linear power spectrum dominates and 
leading non-linear corrections remain negligible, indicating that we can safely rely on lin- 
ear perturbation theory to study inflationary power spectrum. We also briefly comment 
on the infrared and ultraviolet behaviour of the non-linear corrections. 
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1 Introduction 



It is now widely accepted that primordial cosmic inflation in the early universe [T] is the leading 
candidate to provide both the initial conditions for successful hot big bang universe [2], and a 
natural mechanism of generating the primordial perturbations. The quantum fluctuations of 
the scalar field which drives inflation, the "inflaton" field, are stretched to super-horizon scales 
during inflation and become the seeds of temperature anisotropies of the cosmic microwave 
background and large scale inhomogeneities [3]. The computation of the power spectrum of 
these primordial perturbations has become a well established subject [I]. 

Although this linear picture is well studied, its extension to include non-linear effect has 
not been studied seriously. Only recently, computations to second order metric and matter 
perturbations were carried out [5], reporting divergent behaviors of non-linear perturbations 
dominating over linear order ones. This implies that perturbation theory breaks down in the 
seed generation stage. However, this is not likely to be the case considering the fact that the 
observed power spectrum on large scales, which is supposed to reflect the primeval behavior of 
quantum fluctuations, is as small as 10 -9 [6]. 



We should note that to consistently study the effects of non-linear corrections to the power 
spectrum, one needs third order perturbations. This means the second order perturbation of 
previous studies [5] is not sufficient. The non-linear contributions to the power spectrum in 
terms of other correlation functions have been known in terms of the 5N formalism [TJ, but it 
is not clear if those results also exhibit any breakdown of perturbation theory. 

In this note, we carry out explicit computations of the leading non-linear contributions 
to the power spectrum. We consider the simplest but important case of single field slow-roll 
inflation. The metric and scalar field matter perturbations can be described in terms of the 
gauge invariant comoving curvature perturbation TZ. We further take both large scale limit 
and slow-roll approximation, which greatly simplify the calculations yet give the leading order 
result. Our result shows that, contrary to the previous report [5], the linear power spectrum 
dominates completely and the non-linear contributions are negligible. This assures the validity 
of linear cosmological perturbation theory in handling the quantum generation process during 
inflation. 

This note is outlined as follows. In Section [2] we present the non-linear equation of the 
curvature perturbation in the comoving gauge. In Section [3] we solve the equation order by order 
and find the non-linear solutions up to third order. In Section H] we calculate the contributions 
to the power spectrum from non-linear corrections. In Section [5] we present the conclusion. 



2 Equations 



2.1 Setup 

Our starting point is the action of gravity with a minimally coupled scalar field, 




v{4>) 



(i) 



where the action is written using the Arnowitt-Deser-Misner (ADM) metric [8], 



ds 
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Vertical bars denote 3-space covariant derivatives, is the 3-space curvature constructed 
from jij, and the extrinsic curvature 3-tensor is given by 

Kij = 2^( JV< b- + %-7o-,o) • (3) 
Below, we denote the traceless part by an overbar, 

K ij = K ij -- lij K y (4) 

with K = K\ = Y j Kij. 

By varying the action ([1]) with respect to N, N l and 0, we obtain respectively the energy 
and momentum constraint equations, and the equation of motion of as 

if Jfj - -K 2 - i? (3) + \e =0 , (5) 

3 777pi 

K j ilj -lK ]i +—J i =0, (6) 

i_ (n* - N<n%) - Kii* - + ^ + w =0 , (7) 

where the conjugate momentum, energy density and momentum of are respectively given by 

<f>,o ~ N % 



it 



N 



eJ- 
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(n^) 2 + 0i>|,,j (9) 
J* =n^ K . (io) 

To the second and higher order perturbations, in general, we have couplings among the 
scalar, vector and tensor perturbations. In this work, we consider only scalar perturbations 
in a flat Friedmann background model. As the temporal gauge (hypersurface) condition, we 
take the uniform field gauge such that the perturbed scalar field vanishes, 5<p = 0. This is 
the same as the comoving gauge in the single component case which sets Jj = 0. We take a 
spatial gauge condition which fixes the spatial gauge degree of freedom completely [9|: under 
this gauge condition we have 7^ = a 2 (l + 2H)5ij, where 1Z is related to the perturbed part of 
the spatial curvature . Under our temporal comoving gauge, we call 1Z as the "comoving 
curvature perturbation" . These gauge conditions fix the gauge degrees of freedom completely 
even to non-linear orders pU]. To summarize, our gauge conditions are 

50=0, (11) 
li3 =a 2 {l + 2n)5 ij . (12) 

(jSJ) and (J7j) together with the trace part of provide a complete set of equations to have 
a closed form equation for the curvature perturbation 1Z. In our gauge, ([7]) and become 

* 4. (k+ * n *n*\ * dv ri^ 
-Jp + [ K + 7p-Jp N hv = a0' (13) 



- 2 -K 2 +(^ + v) =R^ , (14) 
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where an overdot denotes a time derivative with x° = t. By removing K in these equations we 
can derive a relation between TZ in and N. Further, fl3]) and ([6]) become 

K~(2N i ]i -^) , (15) 

%=^#V (16) 

By removing K in these equations, we can derive another relation between TZ in 7^ and N . 
Combining the two relations between TZ and N we can derive a closed form equation of TZ even 
to non-linear orders in perturbations. 



2.2 Non-linear equations 

Now we consider non-linear perturbation theory of TZ [TO]. We can explicitly combine f[T3"j) . 
(fill) . ( TT5|) and ( IT6|) to write the equation purely using the curvature perturbation TZ to all 
perturbation orders. We will write up to third order, since we are interested in the next-to- 
leading corrections to the power spectrum. Although the full non-linear equation of TZ is very 
lengthy, we can gain more control by taking two approximations. We are interested in the 
behaviour of TZ on large scales. Thus we keep the leading correction terms in the large scale 
limit, which include two spatial derivatives. This corresponds to taking the super- horizon limit. 
Further, in the single field inflation model we consider, the slow-roll approximation is valid with 
a tiny deviation of 0(e), where e = —H/H 2 = ) 2 /(2mpj) is the slow-roll parameter. The 
zeroth order terms in e give the leading contributions and we keep only those terms. 
From ffT3"j) . ( 1X3]) . (fTS]) and f lTB]) . up to third order we have 



7Z + 3H7Z- %1Z=\ 



- ~lZ'%i - 21ZMZ - ^A- 1 (TZ^ATZ) . 
13 

+ —1ZlZ> l lZi + A1Z 2 A1Z - TLA' 1 (W l A1Z) . 

+ ^A- 1 (AKWAIZ + ITZ^TZ^W 1 + A1ZA- l Bi + TZ ,ij A' 1 B 2 j) 



where A = 5 l3 didj and A are the Laplacian and inverse Laplacian operators, and 



(17) 



b 2 = -n^n i + 21ZAK - -a- 1 hz^An) . . (is) 

We can write (|T7|) in the Fourier space by introducing the Fourier component of TZ as 

d 3 k 



R = 1 Wf elkxlZk ■ (19) 
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Then, we can find 



K k + 3HTZ k + ^TZ k =1 J ^^^(k - qu)n qi n q2 Ufi + \q x . ^ + jLfc . ^ 



1 

+ ^7 



d 3 q 1 d 3 q 2 d 3 q 3 



5 {3) (k- qi23 )n qi n q2 n q3 



X 



A 2 13 



1 

+ 2 



Q? k ■ Q23 



123 



k 2 



+ 



k q 1 q 1 - q 23 

k 2 



k ■ Q3 <& 

<?1 ■ <?2 + -^-923 ■ <l2 

T23 



k 2 



?23 
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—q 2 ■ q 3 - 2q 3 + 

4 2g 23 



2-<?23 • <?2 

(20) 



where we have introduced a shorthanded notation q\ 2 ... n = q\ + q 2 + • • • + q n . The non-linear 
parts of (1171) and (1201) are valid to the leading order in the large scale approximation. We note 
that the second and third order terms in the right hand sides of ffTTj) and (120]) have A/a 2 order 
factor, thus suppressed in the large scale limit. 



3 Solutions 

We can find the solution of ( 1201) by perturbative expansion. We first consider the linear solution 

(i) 



1Z k > of (J2DD, which satisfies 



On large scales, we find 



it^ + 3HTZ^+ 5^i 1) = - ( 21 ) 

=Ci + Q, (22) 

with C\ and C 2 being integration constants depending only on k. In an expanding phase 
we are interested in, the constant C\ is (relatively) the growing solution, while C 2 /a 3 is the 
decaying one. Thus, we neglect the transient decaying solution when we feedback the linear 
solution to obtain non-linear contributions. The coefficient C\ will be determined from quantum 
fluctuations as we will see in the next section. 

With the linear solution TZ^ = Ci(fe), we perturbatively expand the full non-linear solution 
lZ k in terms of momentum dependent symmetric kernels a£| 

= n k ] + n k ] + n k ] + ■ ■ ■ , (23) 

where Qi{k) = 1. Using (|23|) . the original non-linear equation (|20|) is reduced to simple differ- 
ential equations of the kernels Q n order by order. 
Plugging (I23I) into (|20|) . at second order we have 

g 2 + mg 2 = 1 ^2g 2 2 + ~ gi . q2 + ^k- Ql \ . (24) 



4 See also Ref. [TT] for the perturbative solution of the density contrast S in Einstein-de Sitter universe. 
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This can be solved to give 

^2(91,92) 



2a 2 H 2 



7 q 2 
H + jQi ■ 92 + ■ Ql 



(25) 



The symmetrized kernel Q% (91, 92) which we will use to find higher order kernels is then given 
by exchanging the arguments, 

1 



^ s) (9i,92) = ^[£2(91, 92) +^2(92,91)] 



With the second order kernel, from 



we can write the second order solution 7c£' as 



K 



(2) 
k 



(26) 



(27) 



For the third order kernel 3 we have 

1 

a 



3 + 3HQ 3 =— 



^2 S) (92,9s) 



2 (ql + qls) + \<h ■ 923 + J^fc ■ q 23 + f|fc ■ q x 



2 

23 
2k 2 ' 



+ 



, > 13 2g? fc • q 3 q 2 

- A % ~ -Tr<l2 ■ 9s - 7T fc • 92 Tn-Qi ■ 92 + —923 • 92 

2 k 2 k 2 q 23 



+ 



1 ( q\ k ■ 923 . *: • 9i 9i • 923 \ / 7 2 3g| 



2 Vfe ^ 2 



+ 



k 2 ql 3 



-^92 ■ 9s - 2g 3 + ^"92 .; • Q2 



23 



Here, let us split the third order kernel Q 3 into two parts, 

£3(9l, 92, 93) = 031 (9l, 92, 93) + £32(91, 92, 93) 



(2f 



(29) 



where Q 3 \ denotes the terms multiplied by the second order kernel Q<f* and Q 3 2 the rest. For 

is) 

^31, we can find that Q 3 \ is a product of two Q 2 s as 



03i (9i, 92, 9s) = -02 S) (92, 9s)0 2 S; (9i, 923) • 



(«)/ 



(30) 



We can also find 32 as 
032(9i,92,9s) = 



2a 2 H 2 



13 



-4?5 - y 92 ■ 93 k2 



2( & h n fc '^„ „ , ql 
k ■ q 2 - 



k ■ q 3 

—75-91 • 92 + ^923 • 92 
K <?23 



+ - 



1 f ql k - 923 , kqxqi- q 23 \ ( 7 ^ 2 3gf 



2 Vfe ^ 2 



+ 



^ 2 ?23 



-92 • 9s - 2g 3 + tti-923 • 92 

4 ^<?23 



Then the symmetrized third order kernel is written as 



03 S) (9i, 92, 93) = 7y [03(9i, 92, 93) + 5 permutations] 



1 



and the third order solution is given by 



K 



(3) 



^y^ 3) (fc " 9123)0^(91, 92, 93)<7«) . 



(31) 



(32) 



(33) 
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4 Power spectrum 



Having found the Fourier mode solution TZ^, we can write the power spectrum as 



(n kl n k2 ) = (n k i y k i i + 



2tt 2 



= (27r) 3 ^(fci + fc 2 )-^^(A; 1 ) 



<27r) 3 5 (3) (fc 1 + fc 2 



2tt2 



,(12) 



(22), 



>( 13 ), 



+ 



(34) 



Here we have grouped the terms of the same perturbation order. (13"4"|) is diagramatically shown 
in Figure[TJ 



• • — = — • • — + —•-{*"•— + —<;;;;>— + -W- * 

T^) ^(fc) ^>(*) 7>£ 3 >(lfc) 



Figure 1: Diagramatic representation of the power spectrum up to next-to-leading order cor- 
rections. Up to this order, all the corrections include one internal momentum loop integral 

f 12) 

so that they can be dubbed "one- loop" corrections. However, note that , which includes 
bispectrum, is of higher order than the other two corrections. 

Before we proceed, we make some remarks. As mentioned in Section [I] the one-loop correc- 
tions to the power spectrum are known from the SN formalism [7] . Especially, the structure of 
these corrections in terms of correlation functions which we show below is precisely the same, 
as it should be. However, there are two differences. First, we obtain our results by solving 
the non-linear Einstein equation. This provides an alternative approach to non-linearity in the 
power spectrum. More importantly, we present the explicit momentum dependences of these 
correction terms, which are hidden in the derivatives of the number of e-folds N in the SN 
formalism. 



4.1 Linear power spectrum 

In the context of inflationary cosmology, the brackets of (134")) are taken with respect to the 
vacuum state of the perturbation operator, in our case 71. Thus, to estimate the linear power 
spectrum of 71, we can expand 7Z k in terms of the creation and annihilation operators of 
quantum harmonic oscillators, namely, 

T^k = ~ ( a kUk + a ] _ k ul} , (35) 

where z = a<p/H [12] and the creation and annihilation operators satisfy the canonical commu- 
tation relations 

[a k ,a\] =(2n) 3 5^(k-q), (36) 
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otherwise zero. Then, the mode function equation of u k is given by 

u'l+(k 2 --)u k = 0, (37) 



z 



where a prime denotes a derivative with respect to the conformal time drj = dt/a. The solution 
which satisfies the Bunch-Davies vacuum boundary condition in the zeroth order slow-roll 
approximation is known to be [13] 



u k = -L (\ + i^f] e* k ^ — ► ^L^*) . (38) 



As k — > we have the asymptotic solution \u k \ oc z, which can be read from ( 13?]) . This is the 
constant amplitude of the growing solution C\{k) we found in (122|) . Note that using the exact 
solution (!38|) . we can recover the same second order kernel (1251) in the large scale limit and 
slow-roll approximation. 

Using (1331) and (I3"51) and comparing with (13~41) . we have 

« } ) = (2 7 r) 3 ^)(fc 1 + fc 2 )^ = (2vr) 3 5( 3 )(fc 1 + k 2 )^-V^(h) , (39) 

so that we can find the linear power spectrum as 

k 3 K| 2 (H\* ^/n 2 



which is the well-known scale invariant power spectrum produced during inflation [3J. In our 
approximation there is no momentum dependence for . 

4.2 Leading corrections to the power spectrum 

Next, we consider the leading correction to the linear power spectrum, the terms inside the first 
square brackets of (1341) . From (1271) . we have 

«)=J ; ^* ,3, (^-'^l 2 )^ , (< ^ l, g2 )(< , n l »nt , 

d 3 gid 3 g 2 (3) (s) 1 

-^^5^(k 2 - q 12 )g ( 2 >( qi ,q 2 )- 

x ((a kl u kl + a ] _ k u* k ^j (a qi u qi + a ] _ qi u* q ^j (a q2 u q2 + al q2 w* 2 J^ . (41) 

What we can see immediately is that, we have the combinations of three creation and/or an- 
nihilation operators, such Thus, we have one remaining creation or annihilation 
operator after using the commutator relations for two of them, which vanishes since it is sand- 
wiched between the vacuum states. Thus, 

(1W2)\ _ /~,(2W1) 



« ( 42 ) 



Hence, 

p£ 2) = , (43) 

and non-vanishing corrections appear only in the next order. 

Before we move to the next-to-leading order corrections, we consider further. What 
we can first note from (|41l) is that it is sourced by the primordial bispectrum, 

(^XMa) = (2vr) 3 5 (3) (fci + k 2 + fc 3 )ffe(fci, fc 2 , fe 3 ) . (44) 

If the distribution of IZ^ is not Gaussian, in general we have a non-vanishing bispectrum. 
Indeed even for single field slow-roll inflation this is the case [Hj. This intrinsic non-Gaussianity 
comes from the contributions around the moment of horizon crossing, which we do not take 
into account by construction [15] . Moreover, although the leading and the next-to-leading order 
corrections all include one internal momentum integral so that they are usually classified as 

(12) 

one-loop corrections, clearly is of higher order. 



4.3 Next-to-leading corrections to the power spectrum 

Now we consider the next-to-leading order corrections to the power spectrum. As we can see 
from (134"]) . there are two contributions. One is given by the quadratic combination of the second 
order solution, and the other by the product of linear and third order solutions. 

4.3.1 V% 2) 

Using the second order solution (1271) . we can find 

«') =/ ^H^ ,3) ( fe i - ftJtf'ta.*) / t^ 3 ^ - fcOtf'fo.*) 

x W^K^S) . (45) 

The physically relevant correlations can be understood in terms of contractions. Since we are 
correlating different perturbations, we have two different ways of contraction as 

(«*) = <(*S?*8>)(*g?*8>)) + ((«')(«))■ (46) 

They correspond to the connected diagram, the third one in Figured) Meanwhile, the remaining 
contractions are within the same perturbations and thus irrelevant: we are interested in the 
correlation between different perturbations. This corresponds to a disconnected diagram, 



Using the above contractions, we obtain 

vf\k) =ft d' qS ^ q , k - ,)#>(-,, - fc + g f£y v *y - 9|) . (48) 

Z7T J q #6 — q\ 
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As is a constant in our approximation, the two linear power spectra inside can be pulled 
out of the integral. We introduce the magnitude of q and the cosine \x between q and k as 
q = rk (0 < r < oo) and k ■ q = k 2 rfi (—1 < ji < 1). After the angular integrations, we have 



V^\k) 



k 



420 V aH 

x 



(ii) 
n 



dr 



r 2 1 1 — rill + r| 



(-41 + 64r + 14r 2 + 84r 3 + 70r 4 + 14r 5 + 14r 6 - 6r 7 - 6r 8 ) 1 1 - 
+ (41 + 64r - 14r 2 + 84r 3 - 70r 4 + 14r 5 - 14r 6 - 6r 7 + 6r 8 ) 1 1 + ? 

We have overall momentum dependence as oc k 4 . 



(49) 



4.3.2 V% 3) 



Next, we move to the rest two terms. We find that 

d 3 qid 3 q2d 3 q 3 (3 



S (3) (k 2 - g 123 )^ s) (9i, 92, <fe) (K%K<gK<gK<g) . (50) 



We can collect relevant correlations as before using contraction. Here, we are correlating a 
single to one 1Z^\ which contains three T^-^'s. Thus all the possible combinations of 
contractions include cross correlations. That is, 



F=T 



<' («><>) > K<'(« >K «)) + (<'(^ ,R «) 



+ 



K^nw))- 



(51) 



Thus, any correlation automatically includes the meaningful one, i.e. the contractions between 
different perturbations. Then, after some computations, we can find 



,(«) 



(')/ 



V { n\k)V y n\q) 
k 3 q 3 



(ii)/ 



(52) 



This can be also analytically integrated with respect to angles and we find 



K 3 \k) 



128 V aH 



k 



?(ii) 

K 



dr 

~l2 



2r (-71 - 97r 2 - 25r 4 + r 6 ) 
+ (l- r 2 ) 2 (-7-2r 2 + r 4 ) log 



1 — r 



1 + r 



k V r^(ii)^ 2 



V- 



K 



16 \aH 

V { K a \k)+V { £ b \k). 



dr 



2r (-30 - 23r 2 + 3r 4 ) + (2 - 5r 2 + 3r 6 ) log 



1 — r 



1 + r 
(53) 



We have different k dependence for the two terms V^ 3a ^ oc k 4 and V^°' oc k 2 . 



,(136) 
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4.4 Numerical integration 

Now we have to integrate the next-to-leading order power spectra in ( 14"9"|) and (|53|) . As the 
linear power spectrum is scale invariant, i.e. 00 we nave ^jj 2 '' oc V^ 3 ^ oc fc 4 and 

Pi 13b) oc k 2 . Since we have worked in the large scale (super- horizon) limit, we cannot integrate 
over the whole range of q but we have to introduce a cutoff in the maximum of q. As our basic 
perturbation equations are valid only in the large scale limit, our analysis is not valid near and 
inside the Hubble horizon scale kjj = aH. Thus we may set g max = k^, which gives 

r ma x = ~r ■ (54) 
k 

Note that as this bound itself is k dependent, the bound introduces additional scale dependence. 
For infrared side we just take a conservative range 

r min = 10- 1000 ^. (55) 

Although we also have logarithmic divergences of the integrals in r — > limit, the infrared 
cutoff does not affect the result appreciably (but see the discussion below). The resulting power 
spectra under such scale dependent bounds are shown in Figure [2J The result apparently shows 
that non-linear contributions are completely negligible compared with the linear contribution in 
the observationally relevant scales. Note that similar conclusions hold for the power spectrum 
of the field fluctuation in the uniform curvature gauge [16] . 

Our result shows that the effect of leading order non-linear terms to the primordial power 
spectrum due to a single minimally coupled scalar field is completely negligible under our as- 
sumptions of the large scale and the slow-roll. Part of the reason can be found in the [k/ (aH)} 2 
suppression factor for the nonlinear terms in (fT7|) and (120]) . The effect of [k/(aH)] 2 suppression 
factor can be either regarded as large scale suppression occurred during the exponential ex- 
pansion, or rapid decaying in time during the same expansion. That is, during the accelerated 
expansion a comoving scale rapidly becomes super-horizon scale. Since our leading order non- 
linear terms already have [k/(aH)] 2 terms, the non-linear contributions are suppressed during 
the accelerated stage as the evolution drives the comoving scales to outside the horizon. 

We make a brief discussion on the integrands of the non-linear corrections. In the limit 
r — > 0, the leading term of each integrand is 1/r and thus is logarithmically divergent. More 
specifically, if we restrict ourselves to a box of size L = 1/ (aH), we have terms with log (kL) [17J. 
Thus, if we push the size of the box to literally infinity, we face a divergence coming from infrared 
region [18]. Indeed, we have checked that as we push the infrared cutoff towards larger and 
larger scales, the contributions of the non-linear corrections increase. This infrared divergence 
may be removed with appropriate manipulations, e.g. boundary effects [19] . 

( 13) 

Conversely, we can find that the integrals in diverge in the large r limit. This should 
be regarded as the breakdown of our approximation near and inside the horizon scale. Thus 
the ultraviolet cutoff we have chosen in (1511) . horizon scale cutoff, is the maximum value of r 
we can have: beyond this value, we are probing sub-horizon regime where our approximation 
is invalid. We may choose more conservative cutoff, for example IQ^aH/k. The net effect of 
pushing r max to a smaller value is to suppress the non-linear corrections, since as mentioned 
above the integrands in P-j^ become bigger at larger r. 
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' TZ 
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^(22) / 

-p(13o), ^(136) / 
72. ' TZ / 

i i . ./.] 



1Q -35 

1Q -30 1Q -25 1Q -20 1Q -15 1Q -10 1Q -5 

k/(aH) 

Figure 2: Total second order power spectrum and the contributions from non-linear corrections 
are shown for scales between the horizon scale to 10 30 factor larger than the horizon scale at, 
say, the end of inflation; this covers the scales relevant to current observations. Note that V^ 3a ^ 
has negative value. We take the absolute value of V^ 3 ^ which is denoted by a dashed line. 



continues the slope till the large-scale limit. V^' + V^ a> + V^°"' nearly coincides with 



(22) 



,(13o) 



,(136) 



V^ 3b \ and V^ 1 ' nearly coincides with V£ 



,(n) 



(total) 



5 Conclusions 

In this note, we have studied the non-linear corrections to the power spectrum of the comoving 
curvature perturbation produced during single field slow-roll inflation. All the scalar perturba- 
tions in the metric and the inflaton field are described in terms of the gauge invariant comoving 
curvature perturbation 1Z. If 1Z to the linear order is Gaussian, we need up to third order 
perturbation to describe the leading non-linear contributions to the power spectrum. Under 
the assumptions of large scale limit and slow-roll approximation, we have solved the equation 
of TZ perturbatively up to third order. Using these solutions, we have computed the power 
spectrum V-ji including the leading non-linear corrections. The resulting power spectrum is, on 
super-horizon scales, dominated by the linear contribution and the non-linear corrections 

are negligibk^]. Our study indicates that we can safely rely on linear cosmological perturbation 
theory to study power spectrum originated from quantum fluctuations. 

5 The non-linear perturbations in the gradient expansion |20j should be very closely related to our calculations 
in this note. We would like to address this point in a separate report. 
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